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ABSTRACT 

The  current  study  is  aimed  at  determining  optimum  gear  ratios  of  mechanical  driven  systems  using  a  three-step  hevel 
helical  gearbox  and  a  chain  drive.  A  n  optimization  problem  with  the  objective function  as  the  minimum  gearbox  height 
was  constructed.  Six  input  parameters  including  the  total  system  ratio,  the face  width  coefficients  of  the  bevel  and  the 
helical  gear  sets,  the  allowable  contact  stress  and  the  output  torque  were  discussed,  and  their  effects  on  the  optimum 
gear  ratios  were  considered  through  the  implementation  of  a  simulation  experiment.  The  findings  revealed,  the 
evaluation  of  the  factor  influences,  moreover,  the  proposed  equations  indicated  that  the  optimum  gear  ratios  can  he 
obtained  with  high  precision. 
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1.  INTRODUCTION 

Mechanical  drive  systems  are  popular  in  different  fields  of  industry.  They  are  used  to  transfer  motion  and  torque  to 
axes.  Additionally,  their  size,  mass  and  cost  depend  to  a  large  extent  on  the  transmission  ratio  of  the  system  drives. 
Therefore,  calculating  the  optimal  transmission  ratio  of  all  levels  of  the  mechanical  drive  system  has  attracted  great 
attention  from  scientists. 

So  far,  the  gear  ratios  of  mechanical  drives  along  with  the  gearboxes  have  been  generally  calculated  for 
different  gearboxes.  They  are  defined  for  helical  gearbox  [1-19],  bevel  gearbox  [1,  3,  20,  21,  22],  and  worm  gearbox 
[23-26].  Furthermore,  they  have  been  identified  for  gearboxes  with  different  steps  such  as  2-step  [1-8],  3-step  [9-15] 
and  4-step  [16-19]  gearboxes.  In  addition,  the  studies  on  the  determination  of  optimal  transmission  ratio  have  also 
been  performed  with  different  objective  functions  such  as  the  minimum  gearbox  cross  section  area  [9,  12]  the 
minimum  gearbox  length  [7,  11,  25],  or  the  minimum  mass  of  gears  [10,  13,  14,  16,  18].  Recently,  several  studies 
have  been  carried  out  on  obtaining  gear  ratios  of  mechanical  driven  systems  using  a  gearbox  and  V-belt  [5,  22,  27,  28, 
29]  or  chain  drive  [30,  31]. 

This  study  focuses  on  calculating  the  optimum  gear  ratios  of  mechanical  drive  systems,  using  a  three-step 
bevel  helical  gearbox  and  a  chain  drive  to  achieve  the  smallest  gearbox  height. 
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2.  METHODOLOGY 


The  height  of  a  mechanical  driven  system  using  a  three  step  bevel  helical  gearbox  and  a  chain  drive  (figure  1)  can 
be  calculated  as: 

h  =  max(de2l,dw22,dw23,d2,)  (1) 

Where,  de2l ,  d w-,-,  and  r/w23  are  the  pitch  diameters  of  the  first,  the  second  and  the  third  steps,  respectively;  d2  is  the  driven 
sprocket  diameter  of  the  chain  drive. 

As  a  result,  the  optimization  problem  can  be  written  as: 

minimize  h  (2) 

With  several  following  constraints: 

1  <  Ui  <  6 

1<  u2  <9 

\<u3<9  (3) 

1  <  uc  <  6 

W  j  j 

Where,  u  t ,  un  and  U-,  are  the  gear  ratios  of  the  first,  the  second,  the  third  steps,  respectively;  uc  and  ut  are  the  gear  ratios 
of  the  chain  drive  and  the  total  gear  ratio  of  the  system. 

The  above  equations  reveal  that  the  requirement  for  solving  the  optimization  problem  is  to  determine  de2l ,  dv22 , 

^w23  ' 
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With  the  chain  drive  unit,  the  pitch  diameter  of  the  driven  sprocket  can  be  expressed  as  [32]: 

d2  =  dl-uc  (4) 

In  which,<://  is  the  pitch  diameter  of  the  drive  sprocket  which  is  determined  by [32]: 

c/j  =  p/sm(n/z1)  (5) 


Where,  Zjis  the  number  of  teeth  of  the  drive  sprocket;  Zican  be  found  by  the  following  regression  equation  [30]: 

z,  =  32.4  -  2.4  -uc  (6) 

P  is  the  chain  pitch  (mm)  determined  based  on  the  value  of  the  design  power  capacity  P  which  is  considered 

as[32]: 

P  =  Prk-kz-kn  (7) 

In  which,  P,  is  the  chain  drive  power  rating  (kW)  which  can  be  determined  by: 

P1  =  7j-nI/(9.55-106)  (8) 

Where  in,  n ,  is  the  drive  sprocket  revolution  (rpm): 

«i  =  nm  /  ug  (9) 


Ti=Tout/{ucrjcrib) 


(10) 


In  which,  Tjc  =0.95=0.97  is  the  chain  drive  efficiency  [32]);  TJh  =0.99-4). 995  is  the  bearing  efficiency  [32]);  7j 
and  Tnut  are  the  drive  and  the  output  torques  (Nmm).  Choosing  tj  =  0.96  ,  7jh  =  0.992  and  substituting  these  values  into 
(10)  gets: 


T\  =  1.0502  /  uc 

k,  k .  and  kn  are  coefficients  which  can  be  found  by  [32]: 
k  =  kd-kp-kc-  kadj  ■  klab  ■  kcon 
kz  =  25/Zl 

K  =  «oi '  ”i 


(11) 


(12) 

(13) 

(14) 


In  equations  (12)  (13)  and  (14),  the  coefficients  kd  ,  kp  ,  kc ,  kadj ,  /t|ub  and  kcon  reflect  the  effect  of  the  shock 

factor,  the  drive  position,  the  center  distance,  the  adjusted  possibility  of  the  center  distance,  the  lubrication  conditions, 
respectively;  «01  is  the  tabulated  drive  sprocket  teeth. 
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2.2  Determining  External  Cone  Distance  of  Bevel  Gear  Unit 

The  external  cone  distance  R  of  the  straight  bevel  gear  unit  is  found  by  [32]: 


(15) 


Where,  kR  =  50  (MpaI/3)  is  the  material  coefficient;  [32];  kbe  =  0.25 . . .  0.3  is  the  face  width  coefficient;  [32];  \<JH { ]  is 
the  allowable  contact  stress  (MPa);  KH/jl  is  the  contact  load  ratio  for  pitting  resistance  determined  by  the  regression 
equation  found  from  the  data  in  [32]: 


Km  =  0.25  •  +  0.2  •  &  + 1 .02 


06) 


In  which,  k  =  khe  ■  M,  /  (2  ~kbe ) . 


Tn  is  the  pinion  torque  which  can  be  calculated  by: 

Tu  =  T0Ut/(ug  ■  uc  ■  r)bg  ■  rig  ■  Vc  ■  rib)  (17) 

In  the  above  equation,  Tout is  the  output  torque  (Nmm);  u  is  the  total  gearbox  ratio 577  =0.95. .  .0.97 is  the  bevel 
gear  efficiency  [32] ; 77, ^  =0.96. . .0.98 is  the  helical  gear  efficiency  [32];  1/  =0.95. . .0.97 is  the  chain  drive  efficiency; 
rjb  =  0.99...  0.995  is  the  rolling  bearing  efficiency  [32].  Choosing  r)hg  =  0.96, rjh  =0.97  ,T]c  =0.96  TJb  =0.992  and 
substituting  these  values  into  equation  (17)  gives: 

T1X  =  1.190  •  T0Ut/{ug  ■  uc)  (18) 

Substituting  kR  =  50  and  (18)  into  (15)  gets: 


Re  =  52.984  •  Ju\  +  1  •  3lrout-  j- - r — 

V  1  yj  1  [(l-kbe)kbeUyUgUc-[oH]2] 

The  outer  pitch  diameter  of  the  pinion  of  the  bevel  gear  set  can  be  calculated  by  [32]: 
dell  =  2  -R'/y/l  +  ut 

And  the  pitch  diameters  of  the  bevel  gear  unit  de21is  determined  by: 
de2i  —  ui  '  deu 

2.3  Determining  the  Center  Distance  of  the  Second  Step 

For  the  second  step,  the  center  distance  a w2  is  determined  by  [32]: 

av/2  —  k,n  ’{U2  0  '  yj'f 12  '  !  ([<7"//  ]  '  U2  '  V ba2  ) 


(19) 


(20) 


(21) 


(22) 
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Where,  kHp  =  1 .02  1 .28  is  the  contact  load  ratio  for  pitting  resistance  [32];  K]  is  the  allowable  contact  stress  (MPa); 

In  practice,  [<7^]  =360. .  ,420(MPa);  km  is  the  coefficient  of  material;  km=  43  [32];  Iff  iU!-<  =  0.3 .. .  0.35  is  the 
coefficient  of  wheel  face  width  [32];7]9is  the  pinion  torque  which  is  calculated  by: 

T12  =  Tout/(u2  •  u3  •  uc  ■  Tj\g  ■  r\c  ■  ri3b)  (23) 

Choosing  jjh  =  0.97  ,  TJc  =  0.96  fjb  =  0.992  and  substituting  them  into  equation  (23)  gets: 

7’i2  =  1.134  ■  Tout/(u2  •  u3  •  uc)  (24) 

Substituting  kup  =  1 . 1 ,  km  =43,  kup  =1.1,  u2-u3=  ug  /  m,  and  (24)  into  (22)  gives: 


46.2882(uz  +  1) 


1'out'ul 


[^h]2“S“2-“c'5,6q2 


In  addition,  for  the  second  helical  gear  step,  dwll  can  be  found  by  [32]: 


(25) 


dw22  =2'aw2'U2/  («2  +  !)  ^26) 


2.4  Determining  the  Center  Distance  of  the  Third  Step 

Similar  to  subsection  2.3,  the  center  distance  of  the  third  step  a  ,3  is  found  by  [32]: 


av/3  Km'{U 3+1)'  Tl3  'kHpl[^GH\  '  ^3  '  yf ba3^ 


1/3 


(27) 


Where,  7j3  is  the  output  torque  which  is  determined  by: 

T13  =  T0Ut/(u3  ■  uc  ■  t]hg  ■  rjc  ■  til ) 

Choosing  rj  =  0.97  ,  Tjc  —  0.96  r/b  =  0.992  as  in  subsection  2.3gives: 

T13  —  1.0912  ■  T0Ut/(u3  ■  uc) 

Substituting  (29),  km  =43  and  kH p  =1.1  (as  in  section  2.3)  into  (27)  gets: 


aw3  =  45.6983  ■  (u3  +  1) 


Tout 


I  l<TH]2'u3'uc'1/,ba3 

As  a  result,  the  pitch  diameter  of  the  third  step  can  be  found  by  [32]: 
dv23=2-aw3-U3/{U3+1) 


(28) 


(29) 


(30) 


(31) 


Experimental  Work 

In  this  study,  to  investigate  the  effect  of  input  parameters  on  optimal  partial  gear  ratios,  a  simulation  experiment 
was  designed  and  implemented.  In  the  experiment,  a  2-level  full  factorial  design  and  6  input  parameters  were  selected 
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(table  1).  Thus,  the  experiment  was  designed  with  26  =64  numbers  of  tests.  In  addition,  to  perform  the  experiment,  a 
computer  program  was  built  based  on  Equations  (2)  and  (3).  The  different  levels  of  the  input  factors  and  the  output 
responses  (the  optimum  gear  ratios  of  the  second  and  the  third  step  ( u 2  and  u3 )  and  the  chain  drive  uc  are  presented  in 
table  2. 


Table  1:  Input  Parameters 


Factor 

Code 

Unit 

Low 

High 

Total  gearbox  ratio 

u, 

- 

40 

70 

Coefficient  of  the  face  width  of  bevel  gear  set 

Kbe 

- 

0.25 

0.3 

Coefficient  of  wheel  face  width  of  step  2 

Xba2 

- 

0.3 

0.35 

Coefficient  of  wheel  face  width  of  step  3 

xba3 

- 

0.35 

0.4 

Allowable  contact  stress 

AS 

MPa 

350 

420 

Output  Torque 

T 

out 

Nm 

102 

104 

3.  RESULTS  AND  DISCUSSIONS 

3.1  Influences  of  Input  Parameters  on  the  Optimum  Gear  Ratio  of  the  Chain  Drive 

To  evaluate  the  influences  of  the  input  parameters  on  the  optimum  gear  ratio  of  the  chain  drive  u  ,  figure  2  describes  the 
graph  of  the  main  effects  for  uc .  It  is  noticeable  from  the  graph  that  ur  is  greatly  influenced  by  the  total  system  gear  ratio. 
In  addition,  it  depends  to  a  lesser  extent  on  y/ha2 ,  ,  and  kbe  in  descending  order.  Also,  it  is  not  determined  by  the 

coefficient  of  the  face  width  of  bevel  gear  set  khe  and  the  allowable  contact  stress  AS. 


Table  2:  Experimental  Plans  and  Output  Response 


Run 

Order 

Center 

Pt 

Blocks 

ut 

Kbe 

Xbal 

Xba2 

AS 

(MPa) 

Tout 

(Nm) 

Uc 

U2 

u3 

1 

1 

1 

40 

0.25 

0.3 

0.35 

350 

10000 

4.20 

4.35 

2.27 

2 

1 

1 

40 

0.3 

0.35 

0.4 

350 

100 

4.00 

4.11 

2.27 

3 

1 

1 

70 

0.25 

0.3 

0.4 

420 

10000 

5.60 

5.60 

2.60 

4 

1 

1 

70 

0.25 

0.3 

0.4 

350 

100 

5.60 

5.27 

2.60 

5 

1 

1 

70 

0.25 

0.35 

0.4 

420 

10000 

5.70 

5.87 

2.46 

6 

1 

1 

40 

0.25 

0.35 

0.4 

420 

10000 

3.90 

4.32 

2.33 

63 

1 

1 

70 

0.3 

0.35 

0.35 

420 

100 

5.80 

5.81 

2.32 

64 

1 

1 

70 

0.3 

0.35 

0.4 

420 

10000 

5.70 

5.87 

2.46 

Main  Effects  Plot  for  uc 
Fitted  Means 


Ut 

KDe 

£5 

Tout 

/ 

/ 

40  70  0.25  0.30  0.30  0.35  0.35  0.40  350  420  100  10000 


Figure  2:  Main  Effects  Plot  for  uc 
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Term 


AD 

C 

AC 

AB 

B 

CD 

BD 

BF 

E 

BE 

DE 

DF 

BC 

CF 

AE 

EF 

F 

AF 

CE 


Pareto  Chart  of  the  Effects 

(response  is  uc,  a  =  0.05) 


Factor 

Name 

A 

ut 

B 

Kbe 

C 

Xba2 

D 

Xba3 

E 

AS 

F 

Tout 

0.0  0.2  0.4 


0.6  0.8  1.0 
Effect 


1.2  1.4  1.6 


1.8 


Figure  3:  Pareto  Chart  for  u 


Coded  Coefficients 


Term 

Effect 

Coef 

SE  Coef 

T-Value  P-Value 

VIE 

Constant 

4.900 

0.000 

*  * 

Ut 

1.6000 

0.8000 

0.0000 

*  * 

1.00 

Kbe 

0.05000 

0.02500 

0.00000 

*  * 

1.00 

Xba2 

0.05000 

0.02500 

0.00000 

*  * 

1.00 

Xba3 

-0.2000 

-0.1000 

0.0000 

*  * 

1.00 

ut'Kbe 

-0.05000 

-0.02500 

0.00000 

*  It 

1.00 

ut*Xba2 

0.05000 

0.02500 

0.00000 

*  * 

1.00 

ut*Xba3 

0.10000 

0.05000 

0.00000 

*  * 

1.00 

Model  Summary 

S  R-sq  R-sq(adj)  R-sq(pred) 

0  100.00%  100.00%  100.00% 

Figure  4:  Estimated  Effects  and  Coefficients  for  uc 


The  Pareto  chart  of  the  standardized  effects  for ;/  is  demonstrated  in  figure  3.  From  this  figure,  it  can  be  learned 

that  the  bars  that  represent  input  parameters  including  the  total  gearbox  ratio  (factor  A),  the  coefficients  of  wheel  face 
width  of  the  third  (factor  D),  the  second  (factor  C)  and  the  first  (factor  B)  gear  units  and  the  interactions  AD,  AC  and  AB 
cross  the  reference  line  (with  value  0  on  the  x-axis).  Consequently,  these  factors  are  statistically  significant  at  the  0.05  level 
with  uc . 

Figure  4  shows  the  estimated  effects  and  coefficients  for  uc.  It  is  found  that  u  ,  khe ,  yshtn ,  Wi,m  anc*  t'le 
interactions  ut  ■  kbe ,  ut  ■  Xf/ba ,  and  ut  ■  l//ba3  are  significant  affecting  factors  on  u  .  Hence,  the  relation  between  uc  and  these 
factors  are  determined  by: 

uc  =  5.8-  0.00000001  -ut  +4.667  kbe  -  2.667  -i/rba2  -11.33- ^,3  -  ^2) 

-0.06667  ■  ut  ■  kbe  +  0.06667  ■ul  ■  y/ba2  +  0.1333  ■  ut  ■  y/ba3 
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The  values  of  adj-R2  and  pred-R2  are  100%  (figure  4),  which  indicates  that  Equation  (32)  matches  the  experiment 

data.. 

3.2  Effect  of  Input  Parameters  on  the  Optimum  Gear  Ratio  of  the  Second  Step 

It  is  interesting  that  the  optimum  gear  ratio  of  the  first  step  of  the  gearbox  does  not  depend  on  the  input  parameters 

and: 

M,  = 1  (33) 

That  means,  with  the  objective  of  the  optimization  problem,  the  bevel  gear  unit  is  only  capable  of  transmitting 
torque  without  executing  the  deceleration  task. 

3.3  Effect  of  Input  Parameters  on  the  Optimum  Gear  Ratio  of  the  Third  Step 

The  effects  of  the  input  parameters  on  the  optimum  gear  ratio  of  the  third  stepw3  are  presented  in  figure  5.  It  can 
be  seen  from  the  figure  that  the  largest  effect  belongs  to  ut ,  followed  by  Xf/btP  ,  y/ha3 ,  Toul  and  AS.  Besides,  kbe  does  not 
effect  on  u  ,  as  the  slope  of  its  graph  is  0. 

The  Normal  Plot  of  the  standardized  effects  is  described  in  figure  6.  From  this  graph,  it  can  be  realized  that  un 
largely  depends  on  the  system  gear  ratio  ut  .  However,  it  slightly  depends  on  If/ h  2 ,  Xf/b  3,  T ut ,  AS  and  the  interactions 
«,  •  Vbai  •  u,  ■  Vba3  and  M,  •  Tou,  in  descending  order. 

Moreover,  ut  (factor  A),  y/ha2  (factor  C),  T  (factor  F)  and  the  interactions  ut  ■  Tout  (AF)  and  ur  ■  1f/hn3  (AD)  have 
a  positive  standardized  effect.  If  these  factors  go  up,  u3  increases.  In  contrast.  If/ hc3  (factor  D),  AS  (factor  E)  and  the 
interaction  ut  ■  Ybai  (AC)  have  a  negative  standardized  effect.  That  means  if  these  factors  rise,  u;,  declines. 

Figure  7  demonstratesthe  estimated  effects  and  coefficients  for  u , .  It  was  reported  that  u t ,  If/ br! , ,  y/ba3,  AS,  Tout 
and  the  interactions  ut  ■  If/ hn 2 ,  ur  •  lf/hll  ,  and  ut  -  T  t  are  the  significant  affecting  factors  on  u2  ■  Therefore,  ,  can  be 
calculated  by: 

u2  =  1.679  +  0.07583  -ut  -  8.833  -kbe  +  0.833  -Ybal  +  5.333 -i//ba3  +  0.03333  -u  ■ kbe  - 

-0.03333  •  u,  -if/ba2  -0.13333-m,  -y/ba3  +20-kbe  -y/ba2  (34) 

Equation  (34)  is  used  to  find  the  gear  ratio  of  the  second  step  u2 .  After  gaining  ux  (Equation  33),  un 
(Equation  34)and  u  (Equation  32),  the  gear  ratio  of  the  third  step  can  be  expressed  by: 

u3  =  ut  /(m1  -u2 -uc)  (35) 

As  tq  =  1 ,  formula  (35)  is  rewritten  as: 


u3=ut/{u2-uc) 

(36) 
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Main  Effects  Plot  for  u2 
Fitted  Means 


Figure  5:  Main  Effects  Plot  for  Optimum  Gear  Ratio  of 


Normal  Plot  of  the  Standardized  Effects 
(response  is  u2,  a  =  0.05) 


Effect  Type 
•  Not  Significant 
■  Significant 


Factor  Name 

B  Kbe 

C  Xba2 

D  Xba3 

E  AS 

F  Tout 


Coded  Coefficients 


Term 

Effect 

Coef 

SE  Coef 

T-Value 

P-Value 

VIF 

Constant 

4.63750 

0.00169 

2751.41 

0.000 

ut 

0.72500 

0.36250 

0.00169 

215.07 

0.000 

1.00 

Kbe 

-0.02500 

-0.01250 

0.00169 

-7.42 

0.000 

1.00 

Xba2 

0.22500 

0.11250 

0.00169 

66.75 

0.000 

1.00 

Xba3 

-0.10000 

-0.05000 

0.00169 

-29.66 

0.000 

1.00 

ut*Kbe 

0.02500 

0.01250 

0.00169 

7.42 

0.000 

1.00 

ut*Xba2 

-0.02500 

-0.01250 

0.00169 

-7.42 

0.000 

1.00 

ut*Xba3 

-0.10000 

-0.05000 

0.00169 

-29.66 

0.000 

1.00 

Kbe*Xba2 

0.02500 

0.01250 

0.00169 

7.42 

0.000 

1.00 

Model  Summary 

S  R-sq 

R-sq(adj) 

R-sq(pred) 

0.0134840  99.90% 

99.88% 

99.86% 

Figure  7:  Estimated  Effects  and  Coefficients  for  u2 


4.  CONCLUSIONS 

This  paper  introduces  a  study  on  determining  the  optimum  gear  ratios  of  mechanical  driven  systems,  using  a  three  step 
bevel  helical  gearbox  and  a  chain  drive.  In  the  problem,  the  minimum  system  height  was  selected  as  the  objective  function, 
and  six  input  factors  including  the  total  system  gear  ratio,  the  coefficients  of  the  face  width  of  three  steps,  the  allowable 
contact  stress  and  the  output  torque  were  inspected.  To  estimate  the  influences  of  these  factors  on  the  optimum  gear  ratios, 
a  simulation  experiment  was  designed  and  conducted.  Considerably,  several  equations  to  calculate  the  optimum  gear  ratios 
were  suggested. 
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